Do Presnel coefficients exist ? 



D. Felbacq\B. Guizal 1 , F. Zolla 2 
1 LASMEA UMR-CNRS 6602 

Complexe des Cezeaux 
63177 Aubiere Cedex, France 
2 Institut Fresnel UMR-CNRS 6133 
Faculte des Sciences de Saint-Jerome 
13397 Marseille Cedex, France 



Abstract 

The starting point of the article is the puzzling fact that one cannot recover 
the Fresnel coefficients by letting tend the width of a slab to infinity. Without 
using the so-called limiting absorption principle, we show by a convenient 
limit analysis that it is possible to define rigorously the field diffracted by a 
semi-infinite periodic medium. 

I. INTRODUCTION 

Over the past decades, the numerical resolution of Maxwell equations has made out- 
standing advances, thanks to the various rigourous methods developed to solve them. One 
may think for instance, in the theory of gratings, of the integral method, the differential 
method or the method of fictitious sources. The tremendous progress of computer capaci- 
ties has been equally important. Generally speaking, it seems that we now have very efficient 
tools to solve any diffraction problem (the 3D is merely a problem of computer power, and 
not a theoretical one). Nevertheless, quoting R. Thorn, the famous mathematician, Field 
medalist and inventor of Catastrophe Theory, "predicting is not explaining". And although 
we are able to compute almost everything, we cannot explain much. This paper points out 
that a simple problem, here a one dimensional one, may lead to deep questions and some- 
what interesting answers. Our starting point is a basic remark concerning diopters. The 
diffraction of a plane wave, in normal incidence, by a diopter separating the vacuum from a 
homogeneous medium of optical index v gives rise to Fresnel coefficients as is well known. 
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This is indeed a fundamental result, commonly used by the opticists. Unfortunately 
diopters do not exist and, consequently, the famous Fresnel coefficients should not exist 
either. The problem may be formulated in the following way: let us consider a slab of index 
v and width h. This gives rise to a reflection coefficient r^. A natural question is: do we 
get r when letting h tend to infinity in r^? When we consider a lossless material {y is 
therefore real), the answer is no, obviously So, how can one measure Fresnel coefficients 
as there is no such thing as a semi-infinite medium in practice ? Two kinds of answers are 
generally put forward. The first answer that seems to come naturally to many is to evoke 
the absorption, and the fact that transparent media do not exist. Let us say that it is a 
poor answer. Indeed, consider for instance the diffraction of a beam, in optical range, by a 
homogeneous layer (for example, silice for which the refractive index is v — 1.47 + ilCU 5 ), 
let's say of a few millimeters (namely the substrate). We can hardly say that there is any 
absorption in that case. But we have here a medium that is almost thousand wavelengths 
wide. Another way out generally put forward is to evoque coherence length. This argument 
is probably correct for the natural light for which coherence length is about micrometer, 
i.e. less than the depth of materials commonly used in optics (See |I| for an exhaustive 
review on this subject). However for a beam whose coherence length is greater than a few 
millimeters one has to find another explanation (it is the case for He-Ne laser beam, for 
which coherence length is about 20 cm). Without using the two above quoted arguments, 
we propose a rigorous answer to this problem, in the more general case of a semi-infinite 
periodic medium. More precisely, we show that it is possible to define the field reflected 
by a semi-infinite periodic medium. We demonstrate that there is no need to use the so- 
called limiting absorption principle nor to use explicitly the notion of coherence length. We 
proceed in two stages. First, we study the diffraction of a plane wave by a finite medium. 
In this case, the diffracted field is characterized by a reflection coefficient. In a second stage, 
we study the limit behavior of the reflection coefficient as the width of the medium tends to 
infinity. All the fields that we consider are z-independent and therefore we may reduce the 
diffraction problem to the basic cases of polarization: E// (electric field parallel to z) and 
H// (magnetic field parallel to z). 



II. REFLECTION BY A FINITE ONE-DIMENSIONAL MEDIUM 

We consider a one dimensional structure made of N identical layers (See [0]). A layer 
is characterized by its relative permittivity e (x), which is assumed to be real and positive. 
For convenience up to the section 4 included, the thickness of a layer is equal to 1. The 
structure is illuminated by a wavepacket of the following form: 




where fi is some measure over the set P of parameters (k, 6): 

P ={(k,9) Gl + x 

Explicit forms of \x will be given in section 4. By means of Fourier analysis, the diffraction 
of such a wavepacket by a finite structure may be reduced to the study of the diffraction of 
a plane wave, which we develop in the following paragraph. 
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We consider a plane wave of wave vector ko illuminating the structure under the incidence 
9 ( k is assumed to belong to the (xOy) plane, we denote /3 = k cos9). The total field 
is described by a single variable function Un(x). When the electric field is parallel to the 
z-axis (E// case) U^(x) exp(ik y sin 6) represents the ^-component of the electric field and 
when the magnetic field is parallel to the z-axis (H// case), it represents the z-component 

1 /2 

of magnetic field . Denoting: (3 (x) = k [e(x) — sin 2 (6 1 )] , the total field Un verifies the 
following equation: 

< x < N : (q-\x)U' N y + q-\x)(3 2 U N = (4) 
and the radiation conditions lead to the following equations 

x < : Un(x) = exp(i/3 x) + (k, 6) exp(— ifiox) 

x > N : U N (x) = t N (k, 9) exp(ip x) ^ ' 

with: q = 1 for E// polarization, q = e for H// polarization. Let xi an d X2 be the solutions 
of equation (f|) verifying 

X i(0) = l, xi(o) = o 
Xa(0) = 0, x' 2 (0) = l 

The fundamental matrix of the system is then 

VM — ( X2 ^ 

and the resolvent matrix is: 

R(x,y) = V(y)V-\x) 

It is the matrix linking the value of one solution at point x to its value at point y. The 
monodromy matrix is finally defined as: 



R(0,1) 



xi(i) 



This matrix characterizes a layer as it allows the writing of the matching between the 
boundary conditions. Indeed, taking into account the propagation conditions (|5|), we derive 
the following relation: 



i(3 (l-r N )J \ifi 



which permits to obtain both values of and t^. When dealing with a wavepacket, we 
write the reflected field under the form: 

Our aim is now to study the limit behavior of as tends to infinity. 
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A. Some properties of matrix T. Connection with Bloch wave theory 



The determinant of matrix T is the value of the wronskian of solutions W\ and W2 at 
point 1. As in that case the wronskian is constant, it is equal to 1. Consequently, T is a 
unimodular matrix and its characteristic polynomial is X 2 — tr (T) X + 1. Therefore the 
eigenvalues of T are real if and only if ~ |tr (X) | > 1. This suggests a splitting of the set of 
parameters: 

G={(M)eP,||fr(T)|>l} 
B = {(M)eP,i|MT)|<l} 

A = {(M)eP,fiMT)i = i} 

When (k, 8) G B, the eigenvalues of T are conjugate complex numbers of modulus 1, 
whereas when (k, 8) G G, T has two real eigenvalues 7 (k, 8) , — ^ where, by conven- 
tion, \"f(k,8)\ < 1. If (k,9) G A, then T has an eigenvalue j(k, 9) of multiplicity 2 with 
either 7 (k, 8) = 1 or 7 (k, 8) = —1. We denote by A the subset of A where T or — T is 
the identity matrix. 

Bloch wave theory is the convenient tool when dealing with propagation equations with 
periodic coefficients. Given a periodic medium with period Y, it consists in searching solu- 
tions of Schrodinger or wave equation under the form u-g (if) e l ^ ' 00 , where u-g (if) is a Y- 

periodic function and k belongs to the so-called first Brillouin zone Y' . For one-dimensional 
media, with period Y = [0, 1[, the theory is quite simple, for in that case Y' = [— tt, +ir[. 
Therefore, Bloch solutions write v $ (x) = u f f ) (x)e l< ^ x , where the so-called Bloch frequency 
(f) belongs to Y', and (x + 1) = (x). Thus we have (x + 1) = e^v^ (x). From the 
definition of T this means that e 1 ^ is an eigenvalue of T, and the above remarks shows that 
tr(T) = 2cos0, which provides us with the dispersion relation of the medium and leads us 
to define <fi a function on B by: 

(k, 8) = arccos Qtr(T)^ (8) 

Obviously, for couples (k, 8) belonging to B, it is possible to define a Bloch frequency (f), 
and therefore there exists propagating waves in the medium: such couples (k, 8) thus define 
a conduction band. If (k, 8) belongs to G then there are only evanescent waves and (k, 8) 
belongs to a forbidden band. 



B. Explicit expression of the reflection coefficient 

We suppose that (k, 8) belongs to P\A. Denoting (v, w) a basis of eigenvectors of T we 
write in the canonical basis of : v = (t>i,t>2) , w = (u>i,u>2) • Eigenvector v (resp. w) is 
associated to eigenvalue 7 (k, 8) (resp. 7 _1 (k,8)). It is of course always possible to choose 
(v,w) such that det(v,w) = 1. After tedious but easy calculations, we get r N in closed 
form from (||): 
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t »( k > 8 )= { \™l 1 ?-?f ( 10 ) 

denoting g(xi, 22) = — — , functions / and g are defined by 

if (M)eG g(k,e) = q(y),f(k,e) = q(w) 

rffko^n I 9 ^ d) = g (v) ' ; ^ = g (w) if |g (v)l < 19 (w)l 

1 ' J U(M) = <?(w),/(M) = <Kv) if |?(w)|<|g(v)| 

Remark: We have \g\ < \f \ and in a conduction band / = g -1 so that we always have 
\g(k,9)\<l. 

Let us denote finally denote: 

z 2 - 1 

H r (Z) = f- 



f -g-y , (12) 

H t = (l-g- 1 f) 2 

We immediately see that the reflection and transmission coefficients are obtained through 
H r and if t by: 

r N (k, 9) = H r (i N ) , t N (k, 6) = H t ( 7 N ) . 



III. ASYMPTOTIC ANALYSIS OF THE REFLECTION COEFFICIENT 

The reflection coefficient defines a sequence of points belonging to the closed unit disc D 
of the complex plane. In order to have a clear understanding of the behavior of {rjv, t^}, we 
interpret this sequence as a discrete dynamical system. As N increases, we want to study 
how the orbits of this system spread over D. 

When (k, 9) belongs to G U (A \ A ) the behavior of the dynamical system is trivial 
as {r^v, t^} admits one cluster point, situated on U = {F, \F\ = J^}. Indeed, in that case 7 
belongs to Recalling that by convention I7I < 1, we see that H r vy N ) tends 

to g (k, 9). As the eigenvectors of T are real, we conclude that g (k, 9) belongs to U. The 
second easy case is for (k,9) belonging to A , indeed {r/v,£jv} is constant and equal to 
whatever N. 

These two cases hand the gaps as well as the edges of the gaps. The case of the conduction 
bands, i.e. for (k, 9) belonging to B is much more complicated and interesting. In the 
following we will skip the mathematical rigor and stress on the physical meaning of the 
results. The interested reader will find a complete and rigorous mathematical discussion in 
§• 

Dealing with a couple (k, 9) belonging to B, the eigenvalues of T now belong to U and 
we may write 7 = e 1 ^ . Obviously, = H r (e lN ^ has no pointwise limit as N tends to 
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infinity. From a geometrical point of view, it is easy to show that the image of U through 
H is a circle V (k, 9) passing through the origin and whose cartesian equation writes: 

(x - $te{z }) 2 + (y- %m{z }) 2 = \z \ 2 , with z = (f 1 + g' 1 )' 1 (13) 

Therefore, {r/v,ijv} describes a set of point on V(k,9). As {r^t/v} does not converge 
pointwise, we turn to another notion of convergence, in some average meaning. An easy 
computation in the case \g(k,9)\ < \f(k,9)\, shows that 

+00 

r N (k,9) =g + g^ 2Nk [fe/" 1 )* - fe/" 1 )*" 1 ] (14) 
k=i 

+00 

t N (k,9) =(l-gr 1 h N J2^ Nk (9f~ 1 ) k (15) 



fc=0 



then the reflected and transmitted fields write 
U*(x, y, t) = J^g (k, 9) exp^") dp: 



+00 „ 

+ E / l U,k 9(k,e)Ugf- 1 ) k -(gf- 1 ) k - 1 
k=i Jp 1 

+00 „ 

UUx,y,t)=Yl / 7 (2fc+1)7V (l-^r 1 ) (gf-yexp 
k=i Jp 



exp i(k.r-a;t) ^ 



i(k-r— uit) 



dpi 



Definition 1: We say that a sequence of functions ip^ (k,0) converges weakly towards 
ipoo (k, 9) if lim f„ ip N (k, 9) f (k, 9) dpi = f B ipoo (k, 9) ip (k, 9) dpi for every <p belonging to 

N >+oo 

^(B,/i). 

We want to pass to the limit N — > +00 in the preceding expressions. What we expect is 
some averaging over the set Y' . Clearly, the limit behavior relies on the properties of pi and 
0. Let us define a convenient class of measures. Denoting C# (Y') the space of continuous 
V '-periodic functions, we put: 

Definition 2: A measure pi is said admissible, if exp (iN<p (k,9)) tends weakly to 0. 

Of course, this looks like an ad hoc property as it allows to get directly the limits of 
interest, but indeed this is a correct way of handling the problem, as it can be shown that 
the measures of interest for our problem, i.e. that of physical significance, will prove to be 
admissible. From the above expression, we conclude that 



U*(x,V,t)^ J g(k : 9)e^ k r -^dp: 
U t N (x,y,t)^0 



We can now conclude by collecting the above results. 

Proposition 1: As N tends to infinity, r N converges weakly towards r^ (k, 9) = g (k, 9), 
tjv tends weakly to 0. 
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IV. DIFFRACTION OF A WAVEPACKET BY A SEMI-INFINITE MEDIUM 



Now, let us choose the incident field as a wavepacket of the form ([|) where /i is of one 
of the following forms: 

I : fi = p(k,6)dk® d6,p G L 1 (P , dk ® d6) 
U .fi = p{d)5 k ®d6,pe L x (] -f , +f [ , dd) 
III : ft = p(k)dk® 5 e ,p E L l (H + ,dk<&5 e ) 

These measures define the most commonly used incident fields. Indeed measures of type I 
correspond to a general wavepacket, measures of type II to a monochromatic beam, and 
measures of type III to a temporal pulse. In order to apply the above results we have the 
following fundamental result: 

Proposition 2: Measures of type I, II, III are admissible. 

We can conclude that, for these measures /i, the diffracted field Ufj (x, y, t) corresponding 
to iV layers converges uniformly towards U^x^y^t) given by: 

Ui(x,y,t) = J^r x (k,8)exp(ik(sin8x + cos0y))exp(-iwt)dii W 

That way, we have obtained a rigorous formulation for the diffraction of a wavepacket by a 
semi-infinite medium. 



V. REFLECTION OF A MONOCHROMATIC BEAM BY A LOSSLESS SLAB OF 

INFINITE THICKNESS 

A very special and interesting case is the case of a simple slab of optical index v: the 
reflection coefficient (r^) is therefore well known : 

1 — exp 2il3h 

1 ' dio eX P 

where r<a is the reflection coefficient for the diopter {r<a = ^ — ^). When the slab is filled 

Po + P 

with a lossy material iV is therefore element of C — {R, HM}, and in that case: 

V(k ,6) eP,r h -^r dio (18) 

On the contrary, when we are dealing with lossless materials, v is real. In that case, 
has an oscillating behavior and does not converge (and a fortiori does not converge to r^ ). 
However, if we consider a limited monochromatic incident beam described by a function 

^■mc(*£) ?/) • 

/k 
^yictx+^y) d(y (lg) 
-feo 

where p(a) G L l {\ — ko,ko[,da) and characterizes the shape of the incident beam. We are 
therefore in the case where the measure /i is of the form II (cf. paragraph V) i.e. : 
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H = p(9)5 k ® d6 (20) 

with p E — 7r/2, 7r/2[, d9). In these conditions the diffracted (reflected) field u h (x,y) 
corresponding to a slab of finite thickness h follows as such : 



/feo 
r h (a)p(a)e^ ax+l3(a ^ da (21) 
-fen 



fka 

-feo 

and the diffracted field Udio(x,y) corresponding to the diopter : 

r-feo 



/BO 
r dio /i(aM«K ( - Q:r+/3(Q)s/) da (22) 
-fen 



-feo 

and we have the fundamental result : 

f feo 



/K 
roo (a)p(a)e^- ax+ ^ da (23) 
-feo 

It is very easy to calculate the reflection coefficient r m in that case. The monodromy matrix 
is: 



(24) 



cos (Ph) \ sin (ph) \ 
-psm(ph) cos(ph) ) 

so that tr(T) = 2 cos (3h, therefore every (k, 9) belongs to BUA and = 
arccos(l/2Tr(T)) = ph. After elementary calculations we find : 

9= pm = r " / = ( ATfl J (25) 

Taking everything into account, we find the weak convergence = rd% and therefore 

u h — > u dio . 



VI. ELECTROMAGNETIC PEROT FABRY 



In numerous experiments in Optics it appears that the light concentrates round particular 
areas of the overall space, namely the rays. In this paragraph, we establish a link between 
the amplitude associated with a ray and the amplitude associated with the electric field. 



A. Simple Perot Fabry 

It is well known that many reflected rays appear when dealing with the reflection coef- 
ficient of a monochromatic beam by a "simple" Perot-Fabry, we are in the same case as in 
the precedent paragraph, for wich the reflected field corresponding to the finite thickness 
slab Uh(x, y)writes as follows: 

/feo 
r h (a)p(a)e l ^ ax+ ^ a ^ da (26) 
-feo 
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^From equation (pf), the reflection coefficient r^can be expressed as a series in the following 
manner: 



r h = r dlo + Y^ 1 -^ l )^ M ( 27 ) 
1=1 

Consequently, u^can also be expressed as a series (the exponential decreasing of r dio allows 
us to reverse the signs J^and J): 

+00 

u h (x, y) = u dio (x, y) + u h>i (x, y) (28) 

i=\ 

with 

^(rL - l)e^ hl P (a)e^ + ^ da (29) 

•Ao 

23q 2/3 

Finally, introducing the two transmission coefficients £i 2 = — and £21 = — -, we 

find an expression of u^i that the opticist can interpret at first glance: 

u h ,i{x, y) = - f'l r^is^e 2 ^^^)^"^^ da,l>l (3Q) 
Uh,a{x,y) = u dio (x,y) 

Each function Uh,i , / G {0, • ■ ■ , +00} is interpreted as the complex amplitude associated 
with a ray labeled by / (cf. figure 1): 

• the term r^-~ 1 t 12 ti2 corresponds to the amplitude associated with a ray /. 

• the term e 2l/3 ^ hl corresponds to a term of phase which expresses the delay of the 
reflected ray / with respect to the first reflected ray (/ = 0) 

• the minus sign before the integral may appear as suspicious. However, for the opticist, 
the interpretation of this sign is easy. The reflection of all rays are of the same nature 
(reflection of a ray from a medium of index v on a medium of index 1) except for the 
first ray (reflection of a ray from a medium of index 1 on a medium of index v) , which 
implies a change of sign in the reflection coefficient. 

Besides, expression (E3) calls for two remarks : 



1. the optical interpretation of functions Uhj is all the more clear as the following condi- 
tions are better fulfilled: a) the supports of the rays I are actually separated, i.e. the 
function p has a narrow support , b) the incident field is sufficiently slanted and c) the 
depth of the slab is sufficiently large with respect to the wavelength. 

2. in the precedent paragraph, we have shown that Uh tends to u d i as h tends to infinity, 
consequently the series in equation (|2lf ) has to tend to zero when h tends to infinity. 
Therefore, each function Uhj behaves as a corrector for the diffracted field (each ray 
except the first one does vanish when the depth of the slab tends to infinity). 
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B. Generalized Perot Fabry 



We consider now a one dimensional structure made of N identical layers. Writing un{x, y) 
the field diffracted by this structure, we have: 



fco 

u N ( x ,y)= I r N {a)p{a)e i{ - ax+ ^ a)y) da 

-ko 



(31) 



Using the expression of in (0) and the same techniques used above, we find : 

+00 

u N {x, y) = Uoojx, y) + ^ u Ni i(x, y) (32) 

i=i 

with 

/fco 
roo (a)p(a)e^- ax+l3 ^ da (33) 
-fc 

and 

/fco 1- -1 

j 2Nl g(k,6) {gf- 1 ) 1 - (gf- 1 ) 1 ' 1 p{a)e^ + ^ da (34) 
-fco L J 

The same comments as in the preceding section can be made. 



C. Speed of convergence : numerical experiments 

Up to now a special attention has been drawn about theoretical aspects of electromagnetic 
diffraction. More precisely, we have proved that, in any case encountered in physics, the 
function U]sr(x,y,t) which represents the diffracted field converges to a function Uoo(x,y,t) 
which one can easily calculate. But it remains to be seen how the function v,N(x,y,t) 
converges to its limit. For instance, in harmonic regime with a pulsatance u, it is of practical 
prime importance to know the number of layers N v from which one can replace, with a given 
precision r], u^{x,y,uS) by u^x^y^uS). From a theoretical point of view, it is very difficult 
to answer such questions even if , in some respects, we outlined an answer in the precedent 
paragraph. We are thus "doomed" to make numerical experiments and to leave the general 
insights of the Theory. Of course, in this paragraph, we do not aspire to the exhaust iveness 
and we only aim at giving some rough estimates. 

In what follows, we are only dealing with a single Perot-Fabry made of silice for which 
the refractive index is v — 1.47 illuminated by a monochromatic gaussian beam solely 

characterized by the wavelength (A = ), the size Ax w and the height x w of the waist 

UJ 

and the mean incidence angle 6 (cf. figure ??). For such a beam the function p(a) introduced 
above writes as follows: 



p(«) = y ex P 



a — ao 
2Ax„ 



exp (-i(a x w - (3 y w )) 



(35) 
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where a = — sin (6>) and y w = —tan (9) x w . In our numerical experiments, we have 
A 

7T X H 

taken the following values : 9 = — , = 200 and — = 200, where H is the "height of 

4 A A 

observation" (cf. figure ??). For these values, we have drawn, in the same figure \uh{H,y)\ , 
\uoo(H, y)\ 2 and \uh,i(H, y)| 2 as functions of y (shifted of ^f) for the following normalized 

h Axu, 
depths (— = 2, 20, 200, 2000, 5000) and for two normalized waists ( — — — = 5, 50) (cf. figures 

A A 
?? to ?? for a waist of 5, and figures 8 to 12 for a waist of 50). Geometrical Optics predicts 

the following locations of the maxima (cf. figure ??): 



h 

A 


2 


20 


200 


2000 


5000 


2/0 
A 


100 


100 


100 


100 


100 


yi 

A 


102.2 


121.9 


319.5 


2295 


5587 



Finally, we have drawn in figure ?? (resp. fig. ??) the total field map 

\u h (x,y) + u inc (x,y)\ for — — - = 5 (resp. for — — - = 50) and for depth — = 200. 

A A A 

In this example, we see that a width of about 1000 wavelengths (for A = 0.5/zm this means 

a 0.5mm width for the substrate) is necessary to obtain the rays described in classical optics. 

In practical experiments, the width of the substrate is usually of the order of one millimeter 

or more, and therefore the Fresnel coefficients are indeed measured. 



VII. CONCLUSION 



What is classical Optics (or coherent Optics)? Since Maxwell, at the end of the nineteenth 
century, the answer seems to be easy: classical optics is the study of the diffraction of an 
electromagnetic field by a body whose size is very large compared with the wavelength 
(mean wavelength). That means that we have to consider classical optics as a limit for 
small wavelengths of electromagnetic optics. Unlike the doxa, we consider that this limit 
is far from being clear. For instance, it is well known that in electromagnetic optics, the 
diffracted field is very sensitive to the polarization of the incident field. On the contrary, 
for small wavelengths, in a lot of applications, the diffracted field is independent of the 
polarization. In many cases, this remarkable property remains mysterious. In this paper, we 
only dealt with one dimensional problem and we demonstrated that the limit analysis (the 
path from the electromagnetic problem to the optical one) is generally ill posed for a plane 
wave (Fortunately, this case does not exist!) but is well posed for the cases encountered in 
physics. Besides, the limit analysis that we have given above does not take into account 
the roughness of the layers. For the sake of simplicity, we did not think fit to describe 
in a realistic way the process of measuring the field: this would have lead us to use a 
spatial convolution process. These remarks are not at all limitations of our study. On 
the contrary, both these phenomena would have improved the convergence of the involved 
sequences. Nevertheless we do think that this Limit Analysis is far more fundamental 
for a good understanding of classical optics than the usual explanations. Take the case of 
spectacles, which are usually intended to improve the vision. Can one be satisfied by such 
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an explanation as that involving absorption? The challenge is now to extend these results to 
the bidimensional or even tridimensional case, for which the very mathematical formulation 
itself is far from being clear. 
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Figures Captions 

Figure 1: Experimental device, position of the waist and height of observation. 
Figure 2: Location of the maxima foreseen by geometrical optics. 

Figure 3: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 5), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = 7r/4), height of observation 
(H = 100), depth of layer h = 2. 

Figure 4: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 5), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
(H = 100), depth of layer h = 20. 

Figure 5: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 5), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
{H = 100), depth of layer h = 200. 

Figure 6: Wavelength (A = 1) , n — 1.47, size of the waist (Ax w = 5), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
(H = 100), depth of layer h = 2000. 

Figure 7: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 5), position of 
the waist (x w = 200, y w = —200), mean incidence angle {9 = vr/4), height of observation 
(H = 100), depth of layer h = 5000. 

Figure 8: Wavelength (A = 1) , n — 1.47, size of the waist (Ax w = 50), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = 7r/4), height of observation 
(H = 100), depth of layer h = 2. 

Figure 9: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 50), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
(H = 100), depth of layer h = 20. 

Figure 10: Wavelength (A = 1) , n — 1.47, size of the waist (Ax w = 50), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
{H = 100), depth of layer h = 200. 

Figure 11: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 50), position of 
the waist (x w = 200, y w = —200), mean incidence angle (9 = vr/4), height of observation 
(H = 100), depth of layer h = 2000. 

Figure 12: Wavelength (A = 1) , n = 1.47, size of the waist (Ax w = 50), position of 
the waist (x w = 200, y w = —200), mean incidence angle {9 = vr/4), height of observation 
(H = 100), depth of layer h = 5000. 

Figure 13: The total field map \uh(x, y) + u inc (x, y)\ 2 with wavelength (A = 1) , n — 1.47, 
size of the waist (Ax w = 50), position of the waist (x w = 200, y w = —200), mean incidence 
angle (9 = vr/4), depth of layer h = 200. 
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